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The interaction between a straight vortex hne in a superconducting film and a soft magnetic 
nanodisk in the magnetic vortex state in the presence of a magnetic field applied parallel to the 
film surfaces is studied theoretically. The superconductor is described by London theory and the 
nanodisk by the Landau-Lifshitz continuum theory of magnetism, using the approximation known as 
the rigid vortex model. Pinning of the vortex line by the nanodisk is found to result, predominantly, 
from the interaction between the vortex line and the changes in the nanodisk magnetization induced 
by the magnetic field of the vortex line and applied field. In the context of the rigid vortex model, 
these changes result from the displacement of the magnetic vortex. This displacement is calculated 
analytically by minimizing the energy, and the pinning potential is obtained. The applied field can 
tune the pinning potential by controlling the displacement of the magnetic vortex. The nanodisk 
magnetization curve is predicted to change in the presence of the vortex line. 

PACS numbers: 74.25. Ha, 74.25. Qt 



INTRODUCTION 



Pinning of vortices in superconducting films by arrays 
of soft nanomagnets placed in the vicinity of the film is a 
subject that has not been much studied in the literature. 
Most experimental 1] and theoretical 05 S S S 0? 
work carried out so far deals with permanent nanomag- 
nets. For both soft and permanent nanomagnets, pin- 
ning results from the action of the inhomogeneous mag- 
netic field created by their magnetization in the super- 
conductor. However, in the case of soft nanomagnets, 
the magnetization depends on the magnetic fields act- 
ing on it, including the field caused by the vortex itself. 
The question is then how this modifies the pinning po- 
tential. Recently, this question was investigated for vor- 
tices interacting with arrays of soft magnetic nanodisks 
using a very simple model in which the nanodisks are 
approximated by point magnetic dipoles with magnetic 
moments free to rotate 8]. This model was suggested 
by the results of Cowburn et. al. 9], which show that 
in small nanodisks the magnetic state is a single domain 
one, with the magnetization of each nanodisk behaving 
like a giant magnetic moment free to rotate. The calcu- 
lations of Ref 8 show that the pinning potential differs 
considerably from that for a permanent dipole, and that 
it can be tuned by a magnetic field applied parallel to the 
film surfaces. Similar properties are expected for the pin- 
ning potential due to soft nonomagnets in general. This 
paper studies the pinning of vortices by nanodisks in an- 
other magnetic state: the magnetic vortex. This state 
is found in larger nanodisks |9l E3 Q ? and is charac- 
terized by the magnetization vector circulating around a 
nucleus of small dimensions compared with the nanodisk 
radius. The motivation to consider this particular state 
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is that it contains the basic physics of nanomagnets in 
general, but its theoretical description can be approxi- 
mated by a simple analytic model. Here, the interaction 
between one straight vortex line and one nanodisk is cal- 
culated using the Landau-Lifschitz continuum theory of 
magnetism 12j for the nanodisk and London theory for 
the superconductor 13]. The magnetic vortex state is de- 
scribed by the rigid vortex model introduced by Usov and 
Peschany [T^, and by Guslienko and collaborators [T5l |. 
These authors obtained the nanodisk equilibrium mag- 
netization using the variational principle to minimize the 
energy. Usov and Peschany 14] proposed an analytic ex- 
pression for the trial magnetization at zero applied field, 
with the magnetic vortex core at the nanodisk center, and 
the core radius as the variational parameter. Guslienko 
and collaborators extended it to finite fields paral- 
lel to the nanodisk faces. They assumed that the only 
effect of the applied field is to displace the magnetic vor- 
tex rigidly, with the core moving away from the nanodisk 
center. The new trial magnetization is the zero field one 
displaced rigidly with the core, and the core displacement 
vector is a variational parameter. The predictions of the 
rigid vortex model are in reasonable agreement with nu- 
merical calculations for small displacements compared to 
the nanodisk radius Qs] . 

This paper assumes the validity of the rigid vortex 
model for the nanodisk placed in the vicinity of the su- 
perconducting film. The total energy of the nanodisk- 
superconductor system is obtained exactly as a function 
of the magnetic vortex displacement vector. The calcula- 
tions carried out in the paper assume small displacements 
compared to the nanodisk radius, and aproximate the 
the total energy by its expansion to second order in the 
displacement. The equilibrium displacement of the mag- 
netic vortex is obtained by minimizing the approximate 
total energy. It is found that, contrary to what happens 
for the isolated nanodisk, the equilibrium displacement 
is a non-linear function of the vortex line and applied 
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fields. The pinning potential for the vortex line is also 
calculated from the approximate total energy. The most 
important contributiont to the pinning potential is found 
to come from the changes in the magnetization caused 
by the displacement of the magnetic vortex. There is 
also a contribution from the interaction between the vor- 
tex line and the magnetization in the magnetic vortex 
core, but it is negligible, due to the small dimensions of 
the core. It is found that the pinning potential can be 
tuned by the applied field. The mechanism is that the 
field controls the magnetic vortex displacement which, 
in turn, modifies the pinning potential. Applications to 
systems of experimental interest are considered by ap- 
propriate choice of the model pararmeters. The results 
suggest that the pinning potential for these systems can 
be estimated by applying the linear response theory of 
elementary magnetism to the nonodidk. The magnitude 
of the applied field necessary to annihilate the magnetic 
vortex is estimated, from which changes in the nanodisk 
magnetization curve caused by the vortex line can be 
aticipated. 

This paper is organized as follows. The theory for the 
nanodisk-superconductor system is developed in Sec. ^ 
First, in Sec. Ill A[ the rigid vortex model for the isolated 
nanodisk is briefiy reviewed. Then, in Sec. Ill Bl the nan- 
odisk in the vicinity of the seperconductor is considered, 
and the main results of the paper are derived. Finally, in 
Sec, mil the application of the model to systems of exper- 
imental interest is considered, and the conclusions of the 
paper are stated. The Appendix gives the mathematical 
expressions needed in the calculations carried out in Sec. 
IIIBI and a brief review of their derivation. 



II. PINNING POTENTIAL 

The superconductor-nanodisk system is shown 
schematically in Fig. ^ The superconducting film is 
planar, of thickness isotropic and has penetration 
depth A. Its surfaces are parallel to the x — y plane. One 
straight vortex line, with vorticity q = ±1, is present 
in the film at a position defined by with respect to 
an origin in the film top surface. The vortex line core 
radius is ^. A thin magnetic nanodisk, of radius R and 
thickness <C is located above the film, with the 




FIG. 1: Schematic view of the superconducting film with one 
vortex line and a nanodisk placed on top. 




(a) (b) 



FIG. 2: Nanodisk magnetization in the magnetic vortex state 
(qm = 1) with core at nanodisk center (a), and displaced by S 
from nanodisk center (b). Also shown in a) is the cylindrical 
coordinate system with origin at nanodisk center used in Eq. 
©. 

faces parallel to the film surfaces, and center located 
above the film at (x = 0, = 0, z = > 0). A magnetic 
field, H, is applied along the x-axis, that is parallel to 
the nanodisk faces and to the film surfaces. 



A. isolated nanodisk 

First, the rigid vortex model for the isolated nanodisk 
is briefly reviewed. The magnetization is written as 
M(r) = Mss{r)^ Mg being the saturation magnetization. 
The nanodisk energy is given by 

Ed = ^ [ d\ I I' 

cx=x^y,z 

4/<i=^/<^''^-«.-H. (1) 

In Eq. (QJ, the first term is the exchange energy, where 
Rq is the exchange length. The second is the magneto- 
static energy, with cr(r) = Mss{y) • n being the density 
of magnetic charges at the surface of the nanodisk, and 
n is the unit vector in the direction normal to the nan- 
odisk surface. The last term is the energy of interaction 
with the external field, with m = J (i^rM(r). It is as- 
sumed that the nanodisk is isotropic, and that there are 
no volume magnetic charges, that is V • M = 0. The 
rigid vortex model is a variational minimization of Eq. 
(PP). The trial magnetization is written in the cylindri- 
cal coordinate system with origin at the nanodisk center 
and z-axis perpendicular to the nanodisk faces shown in 
Fig. 12^1. For H = Usov and Peschany proposed the 
following expression 

s(p) = qmO {a < p< R) 

s(p) = q^^^e±^^^z (0</><a), (2) 



3 



where 6 = z x p. This equation describes the mag- 
netization vector curhng around the nanodisk center ( 
Fig-EJi ), with counterclockwise ( clockwise) rotation for 
Qm = +1 (—I)- The variational parameter a is the radius 
of the magnetic vortex core. The magnetization outside 
the core is parallel to the nanodisk faces. Inside it has a z- 
component, which can be either positive or negative. The 
equilibrium value of a, obtained by Usov and Peschany, 
is 0.7{RqLz)'^^^ What is important for the pur- 
poses of this paper is not the precise value of a, but the 
fact that it is small compared to R. This follows from 
the assumption that Rq^ <^ R. For H 7^ the trial 
mag netization proposed by Guslienko and collaborators 
[13 is that in Eq. Q displaced rigidly by S from the 
nanodisk center, that is s(| p — S |) ( Fig. ^ ). The 
equilibrium value of S is determined by minimizing the 
energy. The vortex core radius is unchanged as long as 
the displaced core remains within the nanodisk, that is 
for 5 < R — a. The simplest treatment of the rigid vortex 
model assumes that the magnetic vortex displacement is 
small, and obtains Ejj to second order m5/R. The result 
is ^ 



Ed{S)=Ed{0)^M^Vd[^ ^ 



(3) 



2xR^ ^'''MsR' 

where Vd = ttR^L^ is the nanodisk volume, and Sx and 
6y are, respectively, the components of S parallel and 
perpendicular to H. In Eq. ^ the term quadratic in 5 
comes from both the exchange and magnetostatic ener- 
gies, whereas the term linear in 5 comes from the inter- 
action with the applied field. The constant x is the nan- 
odisk linear magnetic susceptibility , and is given by0 



1 

X 



Rpt Lz r, 8-R 1 -, 



(4) 



The first term in Eq. is the contribution from the 
exchange energy, and is small because R ^ Rq. The 
second term comes from the magnetostatic energy of the 
magnetic charge density generated at the nanodisk edges 
by the vortex displacement. The density of these charges 
is given by 

_ {z X 6) - p 

The interaction between the magnetic charges in the vor- 
tex core and in the nanodisk edges is neglected. To cal- 
culate the magnetostatic energy to second order in S/R, 
it is sufficient to use a to first order, that is 

(z X ^) • p 



.(1) 



-Q7 



R 



(6) 



The energy of interaction with the applied field also 
comes from cr(^). The magnetic moment generated by 
the displacement of the magnetic vortex is related to cr, 
Eq. ©, by 



m(<5) 



/ 



(fr a p , 



(7) 



where the integral is over the nanodisk edge surface. To 
first order in S/R, the magnetic moment is obtained by 
replacing a by a^^^ in Eq. 0. The result is is 



mW(,5) = -(/™M,Vd 



(z X S) 
R 



(8) 



By symmetry m(^) does contain terms of second order 
in S/R. 

The equilibrium vortex displacement, obtained by min- 
imizing Ed with respect to ^, is 



z X H 



(9) 



In equilibrium, the magnetic moment is 

me, = m(i)(^eg) = X^DH. (10) 

This result shows that x is the nanodisk susceptibility. 
According to Eq. Q, S^q is perpendicular to H, and de- 
pends only on H scaled by Mg/x- As will be shown in 
Sec. Ill BI Mg/x is also the scale for the combined effects 
of the applied field and vortex line field on the magnetic 
vortex. The approximation to second order in 6/R to 
Ed is valid only for xH/Mg < 1. However, it has been 
used outside this range to estimate the vortex annihila- 
tion field, i^an, at which the magnetic vortex is destroyed, 
and the nanodisk switches to the saturated state 15]. 
This estimate assumes that at = Han the magnetic 
vortex reaches the nanodisk edge, that is 5eq = R. Thus, 
according to Eq. Q, Han = ^s/x- this case the nan- 
odisk magnetic moment is the saturation moment, MsVd 
in the direction of H (see Eq. (|1()|) ). This estimate for 
Han is in reasonable agreement with experiments and 
numerical simulations ^, J^] . 



B. nanodisk-superconductor interaction 

Now the nanodisk is assumed to be in the proximity of 
the superconducting film, as shown in Fig.Q] The equilib- 
rium displacement of the magnetic vortex is recalculated 
by minimizing the total energy of the superconducting 
film-nanodisk system, Et{p^-, to second order in 5/ R. 
In the London limit |1^] 

Et{p^,8) = Ed{8) + Ems{S) + Evm{p^,S) , (11) 

where Ems is the energy of interaction of the nanodisk 
with the superconducting film in the absence of vortices, 
and Eym is the interaction energy of the vortex line and 
the nanodisk. The influence of H on the superconductor 
is neglected. This is justifiable as long as H is smaller 
than the lower critical field parallel to the film surfaces. 

The energy Ems results from the interaction between 
the nanodisk and the magnetic field of the screening cur- 
rent generated by it in the superconductor 1^. As shown 
in the Appendix, Ems can be written as an interaction 
between magnetic charges, 

Ems = IJ "^'"^ J d^r' cT{r)a{r')UMs{r;r') . (12) 
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z=0.5X ^z=O.U 



FIG. 3: a) Interaction potential induced by superconducting 
film between two magnetic charges separated by p, and at 
same height z, compared with the Coulomb potential, a) 
XUms vs. p for d = a and two z values. Inset: mag- 
netic charges above superconducting film, b) Ratio Ums to 
Coulomb potential, pUms, vs. p for z = 0.2A and several d. 



Thus, in the presence of the superconducting film, the 
interaction between magnetic charges is modified by the 
addition of Ums to the Coulomb potential ( see Eq. 
In Fig. |31 Ums is shown for two magnetic charges sep- 
arated by p and at the same height and compared 
with the Coulomb potential. For large distances, Ums 
coincides with the Coulomb potential. Large distances 
meaning p ^ A = 2X^ /d for thin films {d <C A) and 
p A for films with d ^ X. For short distances, Ums is 
considerably smaller than the Coulomb potential. The 
effects of Ems on the magnetic vortex are to change 
the vortex core radius, a, and nanodisk susceptibility, 
X- However these modifications are small. The reason is 
that the dominant contributions to a and x come from 
short distances [l3. [T5l|. where Ums is much smaller than 
the Coulomb interaction. For instance, the contribution 
to x~^ from Ems is found to be about one order of mag- 
nitude smaller than that from the Coulomb interaction, 
Eq. Hereafter Ems is neglected. 
The energy Evm is given by p^] 



E^ 



VM 



-y"dVM(r')- b(r'-pj (13) 

= I dV'a(r')*(r'-pJ. (14) 

where b(r) = — V ^{r) is the magnetic field created by 
the vortex line outside the film. There are two contribu- 
tions to EvM' one from the magnetic vortex core, and 
another from the magnetic charge density at the nanodisk 
edges, cr, Eq. Q, denoted by Ey\^. The contribution to 
Evm from the core is negligible, as justified later. In or- 
der to obtain Ey]^ to second order in it is necessary 
to use a to the same order. The first order term is cr^^\ 
Eq. 0. The second order term is, according to Eq. (|5J, 
given by 



r(2) 



R 



(15) 



The contribution from a*^^^ to Evm, denoted Eifll, is 




FIG. 4: Definition of the displacements Sl, St, and angles 0, 
and 7. 



obtained by noting that <j^^^ can be interpreted as result- 
ing from the uniform magnetization 



M^i) = -qM 
Using Eq. it follows that 



z X S 
R 



where 



4^]] = -T/i,mW. B^(-pJ, 



B^(-pJ = ^ f dVb^(r'-pJ 

Jdisk 

= -qBj_{p^)p^, 



(16) 



(17) 



(18) 



is the average over the nanodisk volume of the compo- 
nent of the vortex field perpendicular to the z-direction 
( parallel to the nanodisk faces ), b^. The argument of 
Bx in Eq. (|T^ is —p^ because is defined with respect 
to the origin shown in Fig.Q whereas in the definition of 
the vortex line field the origin is at the vortex line ( see 
Eq. (EH). Thus 

4^1) = qm MsVd g5±(Pv) ^ , (19) 

where St denotes the component of S perpendicular to 
Pv) ( parallel to — z x p^, see Fig. 0|. The con- 
tribution from calculated in the Appendix. The 
result is 



E 



(e2) 
VM 



MsVnqBiip^ 



i?2 



(20) 



where 6l denotes the component of S parallel to 
Bx(— Pv) (I^iS- ©7 and Bi{p^) has dimension of mag- 
netic field, and is given by Eq. (|Xl2j) . 

Now the equilibrium displacement of the magnetic vor- 
tex is calculated neglecting the contribution from the vor- 
tex core to EvM' The total energy is thus 



Et{p^,S)=Ed{0)^ 
qB±{p^)ST - H5y 



2xR^ 
qBi{p^)SLST 



. (21) 
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The components of 5 parallel and perpendicular to H, 
Sx and Sy^ respectively, are related to Sl and St by (see 
Fig. El). 



5y = 



-Sl cos 6 + St sinO 
-Sl sinO — St cos 6 . 



(22) 



Minimizing Ej- with respect to 6, with the vortex Hne 
held fixed at p^, it follows that 



(Pv) 



R 



qB± — H cos 9 + qmqBiH sin 9 



Hsine 



-qmqBi{qB± 



HcosO) 



R 1 



-Bj 



4f • (23) 



Thus 6eq depends only on the vortex line and applied 
fields scaled by Ms/x-> ^^^^ is non-linear in these fields. 
The non-linearity results from the inhomogeneity of the 
vortex field over the nanodisk volume, which is responsi- 
ble for Bi being non-zero. This solution is only valid if 
B1 <1. Otherwise ^gq, Eq. (|^ . does not correspond to 
a minimum of Et{p^^S)^ Eq. The consequences of 

Eq. are discussed in detail in Sec. IIIII 

Now the neglect of the core contribution to Evm^ de- 

(c) 

noted EyM^ justified. Since a is small compared to i?, 

(c) 

^VM estimated by approximating the core by a 

point dipole located at the nanodisk center with the total 
magnetic moment of the core. According to Eq. (|2j, it 
has only the z-component mi^^ = :^Ms7Ta^Lz{2\n2 — 1). 



Thus, Eyj^ 
son why E^^]^ 



bz{6 — + zqz). The basic rea- 

can be neglected is that it is proportional 
to a^, which is already a small quantity. One effect of 
Ey]^ is to modify the vortex core. Since it depends 
on a, its contribution must be added to Ed in order 
to obtain the equilibrium value of a. However, the ef- 
fect is a small one because El> ~ M'^a^/Lz p^], so that 
E^y^lEu - {Lz/aY{hz/Ms), and (hz/Ms) is smah, as 

shown in Sec. IIIII Another effect of Ey^^ is to displace 
the magnetic vortex from the nanodisk center. The dis- 
placement can be estimated as follows. The force exerted 
by the superconducting vortex on the magnetic vortex 



core is I 



VM 



M^VD{a/R)%z/Msi, where £ is the 



typical scale for variations of bz^ namely £ ^ zq for thin 
films {d <^ A) and ^ ~ A for films with d ^ X. This 
force must be balanced by the elastic force M^VdS/xR^^ 
which gives S/R ^ {a/R)'^{R/£)xbz/Ms, and is smaller 
than the vortex displacement caused by Ey]^ at least by 
a factor a/R. 

The pinning potential for the vortex line is defined as 
the equilibrium total energy for the vortex line held fixed 
at p^. That is Up{p^) = ET{p^,Seq). Thus 

xVd 



-[Bi -2qB^H cos 



2{l -BlY 

-h2qmqBiH sin e{qB^ - H cos 6)] . (24) 
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FIG. 5: a) Fields B± and Bi , in units of (/)o/A^, vs. pv 
for = 2.0A. b) Maximum values of B± and Bi in units of 
(/)o/A^ vs. R. The right-hand scale in Gauss corresponds to 
left-hand one for A = 75 nm. Parameters: ^ = 0.1 A, zq — 
0.2A, = O.IA. 



This result shows that the pinning potential can be tuned 
by the applied field. The mechanism is that, according 
to Eq. (|^!^ . H controls the equilibrium displacement of 
the magnetic vortex which, in turn, modifies Up. 

In order to evaluate the consequences of the above re- 
sults it is necessary to attribute values to the model pa- 
rameters. This is done next. 



III. DISCUSSION 

The results of Sec. Ill Bl are applicable to arrays of soft 
magnetic nanodisks placed on top of low-Tc supercon- 
ducting films, provided that the nanodisks are sufficiently 
far apart to neglect dipole-dipole interactions between 
them, and the vortex density is low enough to neglect 
vortex- vortex interactions. Realistic values for the model 
parameters are the following 3, 0, 0, Permal- 
loy nanodisks: 100 nm < R < 500 nm, Lz ~ 10 nm, 
Ms = 0.8 /cG, and Rq ^ 15 nm. Superconducting films 
at zero temperature: A = 75 nm, ^ = 0.1 A, (i ~ 0.2 — 2 A. 
With this value for A, the scale for the vortex fields is 
(/>o/A^ = 3.6 kG. The distance from the nanodisk cen- 
ter to the film surface is chosen as the smallest possible, 
namely zq = 0.2 A = 2^. This value takes into account 
the existence of an insulating layer between the film and 
the nanodisk, of thickness ~ ^, to avoid the proximity 
effect. It is found that for these parameter values the 
vortex fields B± and Bi are small compared to Ms/x- 
As consequences, Bi is negligible in Eq. the mag- 

netic vortex displacement is small compared with and 
nanodisk radius, and its dependence on the vortex field 
is linear. This is discussed in detail next. 

The fields B± and Bi are shown in Fig . O They de- 
pend only on scaled variables, with A as the length scale 
and ^o/A^ as the magnetic field scale (see Eqs. (|A.4|) . 
(|A.5|) ). Their dependencies on pv are shown in Fig. [SJa. 
Both vanish for pv = 0, and have a maximum at pv ^ R^ 
with Bi smaller than B± by a factor ~ 3. The maximum 
values of B± and Bi depend on R as shown as in Fig. Ob. 
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FIG. 6: a) Equilibrium displacements of the magnetic vortex 
{qm = 1) for , if = 0, = 2. OA. b) Maximum values of 
the equilibrium displacements vs. R for H — 0. Parameters: 
^ = O.IA, zo = 0.2A, = O.IA, Ms = O.220o/A^ 



The right-hand scale in Gauss in this figure corresponds 
to the left-hand side one for ^o/A^ = 3.6 kG. Both B± 
and Bi increase non-linear ly with up to d ~ 2. OA. 
For larger they change little, because the vortex line 
field is generated by currents flowing close to the film 
surface. These results indicate that Bi is small. Using 
X ~ 2, which corresponds to ~ 200 nm, Lz = 10 nm ( 
see Eq. Q ), the maximum values of ^i, obtained from 
the Bi data in Fig . [SId, are (^i)mox ~ 5 x 10~^ for 
d = 2.0A and {Bi)max ^ 1-0 x 10-^ for d = 0.2A. Thus 
5^ <C 1, and can be neglected in the denominators of 
Eqs. (|23j) and (|2^ . The terms with Bi in numerator 
of these equations are also small, because they involve 
the combinations BiH^ and BiB±^ which are in general 
smaller than the other terms. 

Results for Sl^cq and Sr^eq at if = are shown in Fig. 
ini In this case, according to Eq. (|23|) . | SL,eq/^T,eq |= ^i- 
Thus, I SL,eq/^T,eq 1 siucc Bi is Small. The curves in 
Fig.inia show this. These plots also show that ^T^eq/R < 
0.15 for d = 2.0A and Sr^eq/R < 0.06 for d = d.2A. In 
Fig.EJb the maximum values of dL^eq and dr^eq are shown 
as a function of R. Both vary little with R. The reason 
is that, according to Eq. (|23|) . the i?-dependence of SL,eq 
and dr^eq is only through the products xB± and x^i- 
These are nearly independent of i?, because there is a 
cancelation between the increase of x with R ( see Eq. 
(^) and the decrease of B_\_ and Bi with R ( Fig. (Sib). 

The above discussion suggest that, for the parame- 
ter values described above, a good approximation to the 
equilibrium displacement is to put Bi equal to zero in Eq. 
(O. In this case Seq^ is identical to that for the isolated 
nanodisk in the applied field Ht = Bx(— Pv) +H. That 
is, Eq. (gni) with Bi = 0, is identical to Eq. Q with 
H replaced by Ht- In this approximation the pinning 
potential, Eq. reduces to 

Up = -^[^i - 2qB^HcosO] . (25) 

This result is, up to a constant, the magnetostatic energy 
of interaction between the magnetic moment induced by 



Ht in the nanodisk, nigg = x^d^t^ and itself. That 
is 

C7, = -lme,-Hr + ^^. (26) 

This simple result is the just the linear response approxi- 
mation from elementary macroscopic magnetism (l^ ap- 
plied to the nanodisk. It depends only on the nanodisk 
susceptibility, x^ and on the macroscopic vortex field 
acting on it, B^(— p^). Results for the pinning poten- 
tial Up{p^)^ based on Eq. (|25|) . are shown in Fig. [7| as 
two-dimensional plots for characteristic values of H. For 
H = 0, Up has circular symmetry, and is the same for 
vortices {q = 1) and anti- vortices {q = —1), with a de- 
generate minimum located in a circle of radius ^ R 
(Fig. [7^). For H 0^ Up has a non-degenerate minimum 
and also a maximum along the direction of H. The plot 
in Fig.[7|3 corresponds to if < B±^ whereas that in Fig. 
[7t is for H > B±. In the latter case the pinning po- 
tential reduces to Up{p^) = qVDxB±{Pv)H cos 6, which 
is identical to that for a nanodisk with permanent uni- 
form magnetization xH. Since B± <C Mg/Xi the range 
of H values for which the vortex line field infiuences the 
pinning potential is small. Thus, except for small if, 
the spatial dependence of the pinning potential is given 
by B±{p^) cos^, and the magnitude is proportional to if. 
This allows for continuous tuning of Up over a wide range. 
It is useful to compare the magnitude of Up to that for the 
pinning potential for a nanodisk with permanent magne- 
tization, equal to the saturation one, U ~ MsVdB±. For 
if < Mg/Xi Up is smaller than /7, since for if < 5^, 
I Up/U I- xB^/Ms < 1 , and for M,/x > H > B^, 
\Up/U XH/Ms < 1. Using Eq. ^ to extrapolate Up 
to the region if > M^/x, it follows that Up ^ U only at 
H ^ Ms/x- 

The linear response approximation is now used to esti- 
mate the effect of the vortex line on the magnetic vortex 
annihilation field, Han- As in Sec. Ill A[ Han is estimated 
as the value of if for which 5eq = R- In the linear re- 
sponse approximation this occurs for Ht = Ms/x- Since 
B^ <C Ms/Xi Han differs little from that predicted for 
isolated magnetic vortex, that is Han ^ ^s/x- Conse- 
quently, the nanodisk magnetization curve is expected to 
differ little from that for the isolated nanodisk. 

In summary then, for the parameter values described 
earlier, the main conclusions are: i) the vortex line field 
causes only a small displacement of the magnetic vortex, 
ii) the pinning potential is the energy of interaction be- 
tween the magnetic moment induced in the nanodisk by 
the vortex line and applied fields and the fields them- 
selves, calculated according to the laws of macroscopic 
linear magnetism, iii) the vortex line has little effect on 
the nanodisk magnetization curve. 

Now the non-linear dependence of 5eq on the vortex 
line field and applied fields, predicted by Eq. is 
considered. The question is whether or not there are 
parameter values for which Bi is sufficiently large that 
the non-linearity is important and what are its conse- 
quences. One limitation on Bi is that Eq. only 
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FIG. 7: Pinning potential for a vortex line by a magnetic 
nanodisk in the vortex state, in units of eoA = (/)o/(167r^A), for 
d = 2.0A, zo = 0.2A, ^ = O.IA, R = 1.5A, = O.IA, Ms = 
O.10o/A^ a) if = 0. b) if = 0.01(/)o/A^ c) H = 0.1(/)o/A^ 



makes sense if S/R < 1. This effectively restricts Bi 
to relatively small values. One reason is that the de- 
nominators in Eq. (|23|) cannot be too small. Another 
reason is that Bi and B± are not independent (see Eqs. 
(|A.8|) . (|A.12|) ). so that if Bi is not sufficiently small, B± 
is large enough to make 6/R > 1. It is found that, in 
general, Bi < B±^ as exemplified in Fig. [SJ To inves- 
tigate this quantitatively, the equilibrium displacement 
of the magnetic vortex for if = is calculated using a 
new set of model parameters. These are chosen in or- 
der to give displacements larger than the ones obtained 
with the parameters mentioned at the beginning of this 
Section. For the superconducting film the parameters 
are the same as before = O.IA, zq = 0.2 A = 2^), 
except for A which is now chosen as A = 50 nm. In 
this case the fields Bj_ and Bi are unchanged in units 
of ^o/A^, but their values in Gauss change because now 
^o/ = 8 kG. The new parameters for the nanodisk 
are chosen as i? = 3. OA, = O.IX, Ms = OAkG. In 
this case x = 3.0 and Mg/x = 0.13 kG. The results for 
the magnetic vortex displacements are shown in Fig. IHt, 
and compared with those obtained from Eq. with 
^1=0. The differences between the curves for Sr^eq 
with Bi and Bi = are due to B^ in the denomi- 
nator of Eq. The results show that the differences 
is small, even when ^T,eg is a considerable fraction of R. 
For SL,eq^ only the curves for Bi ^ are shown, because 
^L,eq vanishes for Bi = 0. The non-zero (^L,eq shown in 
Fig. IHId results from the non-linearity of Eq. with 



the most important contribution coming from the numer- 
ator. Comparing the curves for ^L,eg and ^T,eg in Fig.|St, 
and using | ^L,eql^T,eq 1= B\^ it follows that B\ < 0.3. 
This shows that the validity of Eq. (|23|) requires small 
values of Bi. 

The non-linear effects in the pinning potential Up^ Eq. 
(P^ . for if = come only from Bf in the denominator, 
and are also small according to the above discussion. The 
properties of Up for the new set of parameters are similar 
to those described earlier. 

The results shown in Fig. [SJa also have consequences 
for the magnetic vortex annihilation field, Han- When 
the vortex line displaces the magnetic vortex by a signifi- 
cant fraction of the nanodisk radius it takes only a small 
applied field to annihilate the magnetic vortex by further 
displacing it to the nanodisk edge. In this case if^n dif- 
fers significantly from that for the isolated nanodisk. To 
estimate Han^ Eq. with = is used. In this case, 
as discussed above, the magnetic vortex is annihilated, 
when Ht ~ ^s/x- This estimate for if^n is shown in 
Fig. IHlb. A strong dependence of Han on the vortex line 
position results. The absolute minimum of if^n occurs 
when the vortex line is at the pinning potential mini- 
mum for if = 0, located at = and ^ — i?, since 
in this case the displacement of the magnetic vortex is 
maximum. Thus, when the vortex line is in equilibrium 
with the nanodisk Han is minimum. This result indicates 
that the vortex line in equilibrium with the nanodisk can 
significantly modify the magnetization curve. 

To conclude then, the magnitude of the displacement 
of the magnetic vortex caused by the vortex line depends 
on the particular values of the model parameters. In the 
cases were the displacement is a significant fraction of 
the nanodisk radius, effects of the non-linear relation- 
ship between the displacement and the vortex line field 
are felt, but are small, and the magnetization curve for 
the nanodisk in equilibrium with the vortex line is pre- 
dicted to differ from that for the isolated nanodisk. In 
all cases, the displacement of the magnetic vortex and 




FIG. 8: a) Equilibrium displacements of the magnetic vortex 
{ Qm — 1 ) for H — 0. b) Magnetic vortex annihilation field 
for vortex line at position (xv,2/v) vs. Xv for constant. 
Parameters: ^ O.IA, zq 0.2A, R 3.0\L^ O.IA, Ms 
0.05(/)o/A^ 
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the pinning potential can be estimated by applying the 
linear response theory of elementary magnetism to the 
nanodisk. This suggests that the pinning potential for 
nanomagnets with other geometrical forms for which the 
magnetic vortex state has been reported J^] can be like- 
wise estimated. 
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APPENDIX: LONDON THEORY RESULTS 

Here some results of London theory for the supercon- 
ducting film, obtained in RefsIH [H are reviewed. The 
objective is to write out the mathematical expressions 
needed to carry out the calculations mentioned in Sec. 
iHBl 



1. nanodisk-screening current interaction 

In the arrangement shown in Fig.^ the nanodisk gen- 
erates in the superconducting film a screening current 
which, in turn, creates a magnetic field at the nanodisk. 
The energy of interaction of the nanodisk with this field 
is given by 

Ems = -\j d\m{v) • b,,(r) , (A.l) 

where bsc(r) is the field of the screening current outside 
the film. The field b^c is can be written as 

bsc(r) = -V*sc(r) , 

$,e(r) = y"dV(M(r')- V')C/Ms(r;r'), (A.2) 

where Ums '^^ given by 

UMs{r;r') = J d^fce^k-Cr^-rl) e-fc(-+-') g(k) , 



9{k) = jj—j 



sinhrd 



Vk^ + A- 



(A.3) 



where denotes the component of r perpendicular to 
the z-direction. Assuming that V • M = 0, and integrat- 
ing by parts, Eq. ()A.1|) can be written as Eq. ([T^ . 



2. vortex line field 

The field created by the vortex line outside the film is 
given by 

b(r) = -V$(r) , 

$(r) = -q^ I ^e^^--e-^^Fx(fe) , (A.4) 

where r is the position vector with respect to an origin 
at the vortex line, and 



Fi{k) = 



[{k + r)e^^ + (A: - r)e-^^ - 2k] 



A:r[(A:-r)2e-^^- (A: + r)2e' 



.(A.5) 



Integrating Eq. I|A.4|I over the direction of k, $ can be 
written as 



{kr^)e-^'Fi{k) , (A.6) 

where Jo is the Bessel function of first kind. Thus, the 
vortex field can be written as 

b(r) = bj_{r±,z)rj_+b^{rj_,z)z , 
hi. = ^k\Mkr^)e-'^^F^{k) , 

Using Eq. irO|) and Eq. it follows that 

= p ^S,ik)k^Mkp,)e-''''>Fiik) , 

Sm = ''f}^ 2sinh(feX./2) _ 



kR 



kLz 



3. superconducting vortex- magnetic vortex 
interaction 

Here the contribution from a^^^ to the vortex line - 
nanodisk interaction is calculated. Using Eqs. ([T^ . and 
(|A.4|) . Ey'^ can be written as 



(e2) 



VM- I dz' / de'cj^^\e') 

«o-L^/2 Jo 

(Pk 



(-'^S) / (A.9) 

where \ p' \= R. According to Eq. (|T^ . a*^^^ is given by 

= q^M.,df sin (/?' - 9') cos (/?' - 9') (A.IO) 

where (5' and 9' are, respectively, the angles between 6 
and k, and between p' and k. Integrating over z' ^ 9' and 
it follows that 



5'^ cos 7 sin 7 



(A.ll) 
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where and 7 is the angle between S and (Fig. Using 

Sl = Scosj and St = (^sin7, Eq. (|A.11() is identical to 



_ 2J2(fei?) 2sinh(feL,/2) 



Eq. poll. 
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